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Abstract:  "Locality"  is  a  fraught  word,  even  within  the  restricted  context  of  Bell's  theorem.  As 
one  of  us  has  argued  elsewhere,  that  is  partly  because  Bell  himself  used  the  word  with  different 
meanings  at  different  stages  in  his  career.  The  original,  weaker,  meaning  for  locality  was  in  his  1964 
theorem:  that  the  choice  of  setting  by  one  party  could  never  affect  the  outcome  of  a  measurement 
performed  by  a  distant  second  party.  The  epitome  of  a  quantum  theory  violating  this  weak  notion 
of  locality  (and  hence  exhibiting  a  strong  form  of  nonlocality)  is  Bohmian  mechanics.  Recently,  a 
new  approach  to  quantum  mechanics,  inspired  by  Bohmian  mechanics,  has  been  proposed:  Many 
Interacting  Worlds.  While  it  is  conceptually  clear  how  the  interaction  between  worlds  can  enable  this 
strong  nonlocality,  technical  problems  in  the  theory  have  thus  far  prevented  a  proof  by  simulation. 
Here  we  report  significant  progress  in  tackling  one  of  the  most  basic  difficulties  that  needs  to  be 
overcome:  correctly  modelling  wavefunctions  with  nodes. 
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1.  Introduction 

This  paper  is  based  loosely  on  the  talk  given  by  its  final  author  at  the  2017  Symposium  on 
Emergent  Quantum  Mechanics  (EmQM17).  That  talk  took  as  its  principal  inspiration  the  following 
questions  (from  the  list  provided  by  the  organisers  to  give  focus  to  the  meeting): 

•  Is  the  universe  local  or  nonlocal? 

•  What  is  the  future  of  scientific  explanation?  Is  scientific  metaphysics,  e.g.,  the  notions  of  reality, 
causality,  or  physical  influence,  obsolete  in  mathematical  accounts  of  the  quantum  world? 

However,  in  view  of  the  conference  also  being  the  David  Bohm's  Centennial  Symposium,  the 
speaker  also  briefly  discussed  some  work  relating  to  this  question: 

•  What  is  David  Bohm's  legacy  for  the  future  of  quantum  physics? 

This  paper  also  concerns  these  two  broad  topics,  but  in  the  opposite  ratio  to  that  in  the  delivered 
talk.  Section  2  of  this  paper  addresses  just  the  question  "Is  the  universe  local  or  nonlocal?"  by 
examining  the  historical  meaning  of  the  term  "local".  (The  second  question  listed  above,  which  was 
addressed  in  the  talk,  and  in  References  [1-3],  is  not  addressed  here).  That  section  expresses  the 
views  of  the  invited  contributor,  not  necessarily  those  of  all  the  authors.  Sections  3-6  concern  Many 
Interacting  Worlds  (MIW),  an  approach  to  quantum  mechanics  inspired  by  Bohm's  original  hidden 
variables  theory  [4,5],  and  which  thus  we  hope  will  be  part  of  the  answer  to  the  last  question  above. 

Having  introduced  the  MIW  approach  in  Section  3,  we  turn  in  Section  4  to  how  it  can  give  rise 
to  nonlocality;  that  is,  how  it  can  offer  a  positive  answer  to  this  (here  abridged)  focus  question: 
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•  Are  nonlocal  connections — e.g.,  "action-at-a-distance" — fundamental  elements  in  a  radically  new 

conception  of  reality? 

However,  as  Section  4  notes,  there  are  still  major  difficulties  faced  in  demonstrating  these 
connections  in  MIW  simulations.  Section  5  introduces  a  new  MIW  technique — a  higher-order 
interaction  potential — to  address  one  of  these  difficulties,  namely,  wavefunction  nodes.  Section  6 
applies  this  idea  to  progress  the  MIW  simulation  of  stable  excited  state  nodes  of  a  quantum  particle 
in  one  dimension.  An  alternative  method  to  address  the  node  problem  is  presented  and  simulated  in 
Appendix  A.  Section  7  concludes  with  a  discussion  of  open  challenges  for  MIW  theorising. 

2.  Is  the  Universe  Local  or  Nonlocal? 

The  answer  to  this  EmQM17  focus  question  hinges,  of  course,  on  what  one  means  by  the  term 
"local"  (assuming  that  "nonlocal"  is  simply  its  complement).  Perhaps  surprisingly,  it  seems  [2,6]  that 
the  word  was  not  used  in  the  context  of  interpreting  EPR  quantum  correlations  prior  to  the  1964  paper 
of  Bell  [7].  In  that  paper.  Bell  proved  his  1964  Bell's  theorem;  to  quote  [7]: 

In  a  theory  in  which  parameters  are  added  to  quantum  mechanics  to  determine  the  results 
of  individual  measurements,  without  changing  the  statistical  predictions,  there  must  be  a 
mechanism  whereby  the  setting  of  one  measuring  device  can  influence  the  reading  of 
another  instrument,  however  remote. 

In  other  words,  some  quantum  phenomena  are  incompatible  with  the  joint  assumption  of 
predetermination  (or  causality;  Bell  used  both  terms)  and  locality  (or  separability;  Bell  used  both 
terms).  In  the  above  quote,  it  is  the  negation  of  locality  that  is  characterised,  in  a  way  consistent  with 
Bell's  earlier  definition  of  locality;  to  quote  [7]: 

It  is  the  requirement  of  locality,  or  more  precisely  that  the  result  of  a  measurement  on  one 
system  be  unaffected  by  operations  on  a  distant  system  with  which  it  has  interacted  in  the 
past,  that  creates  the  essential  difficulty. 

Thus,  Bell  intended  to  be  (somewhat)  precise  about  what  he  meant  by  locality.  Unfortunately  he 
did  not  give  a  general  mathematical  definition,  nor  did  he  define  terms  like  "unaffected"  or  (in  the 
first  quote  and  elsewhere)  "influence". 

In  the  theorem  he  proves  that  the  role  "locality"  plays  is  the  following.  The  assumption  of 
predetermination  means  that  an  arbitrarily  long  time  before  any  measurements  are  performed,  there 
existed  in  the  world  a  collection  of  hidden  variables.  A,  that,  together  with  the  future  measurement 
settings,  determines  the  future  outcomes.  Adding  the  "vital  assumption  [of  locality]"  implies  that  the 
outcome  A  of  one  party  (Alice,  say)  cannot  depend  on  the  settings  b  of  a  distant  party  (Bob,  say),  but 
only  by  her  own  setting  a.  That  is,  in  a  theory  9  with  local  predetermination  of  outcomes,  there  exists  a 
function,  Ag(a,  A)  such  that  A  —  Ag(a,  A),  and  likewise  for  Bob.  From  this.  Bell  was  able  to  derive  his 
famous  theorem,  that  there  exist  sets  of  measurement  on  entangled  quantum  systems  whose  results 
cannot  be  explained  by  any  such  model.  (Note  that  Bell  made  an  implicit  assumption  in  1964,  related 
to  freedom  of  choice,  that  P(A\a,  b)  —  P( A).  Here  we  follow  Bell,  as  such  an  assumption  is  certainly 
necessary  [8].) 

In  the  above  formulation,  "locality"  has  a  precise  meaning  only  in  the  context  that 
predetermination  has  already  been  assumed.  This  is  fine  for  the  purpose  to  which  Bell  puts  it. 
If  one  wanted  to  broaden  the  definition  so  that  it  applied  independently  of  the  assumption  of 
predetermination  then,  it  has  been  argued  [2,6,9,10],  the  natural  reading  of  Bell's  verbal  definition 
would  be  as  follows.  An  arbitrary  theory  contains  initial  variables  A,  which  may  or  may  not  be 
hidden,  and  which  may  or  may  not  be  sufficient  to  determine  all  outcomes,  and  is  described  (in  the 
limited  context  we  are  considering)  by  theoretical  probabilities  Pg(A,B\a,  b,  A).  The  theory  is  local  if 
there  exists  a  function  Pg  ( A  \ a,  A)  such  that 


Vfe,  Pg(A\a,b,A)=Pe{A\a,A). 


(1) 
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and  likewise  for  Bob. 

Others  have  argued  [11-14],  to  the  contrary,  that  this  definition  does  not  fit  with  the  other  use 
Bell  makes  of  the  concept  of  locality  in  his  1964  paper,  which  is  the  first  paragraph  of  his  Section  II. 
There  Bell  implies  that,  according  to  the  EPR  argument  [15],  locality  plus  perfect  correlations  of 
outcomes  implies  predetermination.  This  shows,  if  nothing  else,  that  Bell  did  think,  in  1964, 
that  locality  was  an  assumption  that  had  meaning  prior  to  the  assumption  of  predetermination. 
However,  the  obvious  meaning.  Equation  (1),  does  not  work  in  the  EPR  argument.  Even  one  of 
Bell's  most  ardent  admirers  [13]  was  forced  to  admit  this  [14]: 

It  is  simply  not  clear  how  to  translate  Bell's  words  here  (about  locality)  into  a  sharp 
mathematical  statement  in  terms  of  which  the  EPR  argument  might  be  rigorously 
rehearsed.  . . .  [I]t  must  be  admitted  that  Bell's  recapitulation  of  the  EPR  argument  in  this 
paragraph  leaves  something  to  be  desired. 

Regardless,  is  there  really  a  problem  here?  Bell  does  not  say  that  he  believes  that  locality  plus 
perfect  correlations  implies  predetermination.  He  merely  says,  at  the  beginning  of  that  paragraph, 
"...the  EPR  argument  is  as  follows"  (my  emphasis).  In  the  preceding  (opening)  paragraph  he  is 
even  weaker:  "The  paradox  of  Einstein,  Podolsky  and  Rosen  was  advanced  as  an  argument ..."  (my 
emphasis  again).  In  a  follow  up  paper  in  1971  [16],  Bell  is  weaker  still.  He  gives  three  motivations 
for  the  assumption  of  predetermination  (in  a  section  entitled,  unambiguously,  "Motivations").  An 
EPR-style  argument  is  the  third  motivation;  he  does  not  claim  any  logical  deduction  from  locality  but 
merely  appeals  to  the  intuition  of  the  reader  for  the  reasonableness  of  predetermination  by  hidden 
variables  (see  Reference  [6]  for  details).  The  1969  paper  of  CHSH  [17]  seems  actively  skeptical  of  the 
EPR  argument,  saying  it  "led  them  to  infer  that  quantum  mechanics  is  not  a  complete  theory"  (my 
emphasis).  It  is  thus  clear  that,  at  this  time,  the  EPR  argument  was  far  from  being  regarded  as  a 
rigorous  proof  for  the  necessity  of  hidden  variables.  (For  a  discussion  of  the  extent  to  which  the  EPR 
argument  is  such  a  proof,  see  Reference  [18].) 

The  situation  with  regard  to  the  EPR  argument  changed  rapidly,  at  least  in  Bell's  mind,  after 
he  formulated,  in  1976  [19],  a  concept  that  does  allow  one  to  infer  predetermination  from  perfect 
correlations.  This  was  the  concept  of  "local  causality",  stated  most  succinctly  in  a  later  paper  [20] 

A  consequence  ...  of  "local  causality"  [is]  the  outcomes  [in  the  two  labs]  having  no 
dependence  on  one  another  nor  on  the  settings  of  the  remote  [measurement],  but  only 
on  the  local  [measurement  settings]  and  on  the  past  causes. 

In  the  situation  considered  above,  a  theory  is  "locally  causal"  only  if  there  exists  a  function 
Pg ( A \a,  A)  such  that 

Vfc,B,  Pg(A\a,B,b,\)  =  Pg(A\a,  A),  (2) 

and  likewise  for  Bob.  Note  the  appearance  of  B  as  a  conditional  variable  on  the  left  hand  side,  which 
distinguishes  this  concept  from  "locality"  in  Equation  (1).  Moreover,  this  assumption  obviates  the 
need  to  consider  determinism  at  all — it  leads  directly  to  the  Bell  inequalities  that  quantum  mechanics 
violates.  Thus,  Bell  gave  what  I  have  called  [2]  Bell's  second  Bell's  theorem,  in  1976  [19]: 

Quantum  mechanics. . .  gives  certain  correlations  which  . . .  cannot  be  [reproduced  by]  a 
locally  causal  theory. 

Bell  clearly  (and,  I  think  [2,3],  rightly)  thought  local  causality  to  be  a  more  natural  concept  than 
locality  as  per  Equation  (1),  as  he  never  used  the  latter  concept  again.  Regrettably,  however,  he  did 
not  abandon  the  word  "locality".  Rather,  beginning  even  in  1976  [19],  he  sometimes  used  "local"  as 
short-hand  for  "locally  causal",  and,  a  few  years  later,  was  apparently  convinced  that  local  causality 
was  the  concept  that  he  (and  EPR)  had  always  used  [21]  (for  details,  see  [2]).  However,  at  least  in  his 
final  word  on  the  subject  [20],  Bell  showed  his  preference  unequivocally  for  the  terminology  "local 
causality"  over  "locality". 
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Thus  we  may  return  to  the  primary  question  posed  above — Is  the  universe  local  or  nonlocal?  If  by 
"locality"  one  means  "locally  causal",  the  concept  Bell  promoted  for  most  of  his  career  in  quantum 
foundations  [22],  then  the  answer  (barring  more  exotic  possibilities  such  as  "superdeterminism"  [20], 
retrocausality  [23],  and  the  subjectivity  of  macroreality  [24])  is  that  the  universe  is  nonlocal;  it  violates 
local  causality.  To  avoid  confusion,  we  might  agree  to  say  that  the  universe  is  Bell-nonlocal  [2].  If,  on 
the  other  hand,  one  adopts  the  definition  of  "local"  indicated  by  Bell's  1964  paper  and  commonly 
used  in  text  books  [25,26],  then  the  answer  is  that  we  cannot  say  whether  the  universe  is  local  or 
nonlocal.  Operational  quantum  mechanics  satisfies  this  weaker  sense  of  locality,  simply  because  it 
does  not  feature  signalling  faster  than  light,  and  denies  the  need  for  giving  any  account  for  quantum 
correlations  beyond  an  operational  one.  We  can  only  say  that  the  universe  is  nonlocal,  in  this  strict 
sense,  if  we  make  some  other  assumptions  about  its  nature,  such  as  determinism. 

3.  David  Bohm's  Legacy:  Permission  to  Theorise  Radically  New  Conceptions  of  Reality 

The  most  famous  quantum  theory  which  does  make  the  assumption  of  determinism  is  of  course 
David  Bohm's  [4,5].  Indeed,  Bohm's  theory  was  an  inspiration  to  Bell  who  summarised  the  result  of 
his  1964  paper  in  the  introduction-cum-abstract  as 

[A]  hidden  variable  interpretation  of  elementary  quantum  theory  [4,5]  has  been  explicitly 
constructed.  That  particular  interpretation  has  ...  a  grossly  nonlocal  structure.  This 
is  characteristic,  according  to  the  result  to  be  proved  here,  of  any  such  theory  which 
reproduces  exactly  the  quantum  mechanical  predictions. 

Why  Bell  considered  Bohm's  interpretation  to  be  "grossly  nonlocal",  rather  than  nonlocal 
simpliciter,  is  unclear.  Perhaps  it  was  because  the  theory  is  nonlocal  even  in  situations  where 
there  is  an  obvious  local  hidden  variable  theory,  as  in  the  EPR-correlations  [15],  or  the  EPR-Bohm 
correlations  [27]. 

Unlike  operational  quantum  mechanics,  Bohm's  theory  is  a  precise  and  universal  physical 
theory.  Restricting  to  the  case  of  interacting  nonrelativistic  scalar  particles  for  simplicity  of  discussion, 
it  takes  the  universe  to  be  described  by  a  universal  wavefunction  Y(q),  obeying  Schrodinger's 
equation,  where  q  is  the  vectorised  list  of  the  coordinates  of  all  the  particles.  However,  it  also 
postulates  a  single  point  in  configuration  space,  x,  which  encodes  the  real  positions  of  all  these 
particles.  This  "marvellous  point"  [28]  or  "world-particle"  [29]  has  a  deterministic  equation  of  motion 
x  =  vy(x)  guided  vicinally  by  Y(q).  (Note  that  "vicinal"  is  a  synonym  of  "local"  in  the  latter's 
quotidian  sense,  introduced  here  to  avoid  any  possible  confusion  with  "local"  in  the  technical  sense 
defined  in  Section  2.)  Here,  "vicinal  guiding"  means  that  the  world-particle's  velocity  Vf  (x)  depends 
onT(q)  and  finitely  many  derivatives,  evaluated  at  q  =  x.  However,  vicinal  in  configuration  space 
is  not  vicinal  in  3D  space — the  positions  of  Bohmian  particles  in  one  region  of  3D  space  can  affect 
the  motion  of  an  arbitrarily  distant  particle  if  entanglement  is  present.  Since  it  is  the  position  of 
Bohmian  particles  that  encodes  what  an  experimenter  decides  to  measure,  this  gives  rise  to  the  'gross' 
nonlocality  of  Bohmian  mechanics  which  Bell  noted  in  1964. 

Bohm's  original  proposal  [4,5]  actually  used  a  second-order  dynamical  equation  x  =  aVK-jna]  (x)  + 
a Y (x).  The  3D-vicinal  acceleration  avicinal(x)  is  given  by  Newton's  laws,  involving  inter-particle 
potentials  which  drop  off  with  3D-distance.  The  nonlocal  effects  in  Bohm's  theory  arise  from  a 
separate,  3D-nonvicinal,  quantum  acceleration  (x).  Bohm's  publication  of  an  explicitly  nonlocal 
theory  seems  to  have  made  it  acceptable  for  other  physicists  to  publish  realist  approaches  to  quantum 
mechanics  in  direct  opposition  to  the  Copenhagen  interpretation.  This  included:  de  Broglie  in 
1956  [30],  reviving  his  unpublished  idea  from  30  years  earlier  which  prefigured  much  of  Bohm's 
work  and  used  the  first-order  dynamics  described  earlier;  and  Everett  in  1957  [31],  introducing  the 
relative  state  interpretation,  more  popularly  known  as  the  many  worlds  interpretation. 

Taking  inspiration  from  both  Bohm  and  Everett,  two  of  us  plus  Deckert  introduced  in  2014  [29] 
what  we  called  the  Many  Interacting  Worlds  (MIW)  approach  to  quantum  mechanics.  We  suggested 
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that  it  might  be  possible  to  reproduce  quantum  phenomena  without  a  universal  wavefunction  Y(q) 
(except  to  define  initial  conditions).  In  its  place  we  postulated  an  enormous,  but  countable,  ensemble 
X  =  { Xj  :  j }  of  points  x,  in  configuration  space  (similar  ideas  have  been  proposed  earlier  by  a 
number  of  authors  [32-35],  but  they  considered  a  continuum  of  worlds,  which,  in  our  view,  leads 
to  some  of  the  same  conceptual  issues  that  Everett's  interpretation  faces — see  also  [36]).  Each  point  is 
a  world-particle,  just  as  Bohmian  mechanics  postulates,  and  the  dynamics  is  intended  to  reproduce  a 
deterministic  Bohmian  trajectory  for  each  world-particle.  However,  the  trajectory  of  a  world-particle 
is  guided  not  by  a  wavefunction,  but  by  the  locations  of  nearby  world-particles  in  configuration 
space:  3 cy  —  avicinai(x)  +  avnw(Ny),  where  Ny  c  X  is  the  set  of  world-particles  in  the  vicinity 
(somehow  defined)  of  world-particle  xy,  and  aMiw  is  some  fixed  (/-independent)  function  describing 
the  interaction  of  many  worlds. 

In  the  MIW  approach,  probabilities  arise  only  because  observers  are  ignorant  of  which  world  Xy 
they  actually  occupy,  and  so  assign  an  equal  weighting  to  all  worlds  compatible  with  the  macroscopic 
state  of  affairs  they  perceive,  in  accordance  with  Laplace's  principle.  In  a  typical  experiment,  where 
the  outcome  is  indeterminate  in  operational  quantum  mechanics,  the  final  configurations  of  the 
worlds  in  the  MIW  approach  can  be  grouped  into  different  subsets,  x0  C  X  of  world-particles,  still 
extremely  large  in  number,  based  on  shared  macroscopic  properties  corresponding  to  the  different 
possible  outcomes  o.  In  Everett's  approaches,  these  groups  correspond  to  branches  of  the  universal 
wavefunction,  and  one  has  to  argue  that  the  square  modulus  of  the  coefficient  of  each  branch 
somehow  manifests  as  the  correct  probability  for  the  experimenter.  In  the  MIW  approach,  the 
operational  quantum  probabilities  will  be  equal  to  the  number  of  worlds  in  each  group  divided  by 
the  number  of  worlds  at  the  start. 

4.  Nonlocality  in  the  Many  Interacting  Worlds  Approach 

4.1.  General  Considerations 

Just  as  Bohmian  dynamics  is  vicinal  in  configuration  space  but  gives  rise  to  nonlocality  in  3D 
space,  so  too  is  the  type  of  MIW  dynamics  just  described.  Indeed,  it  must  be  nonlocal  because,  like 
Bohmian  mechanics,  it  is  deterministic.  However,  we  can  be  more  specific  about  how  this  nonlocality, 
in  Bell's  1964  sense,  arises  in  the  MIW  approach. 

Consider  a  branch  Ye  (in  the  Everettian  sense)  of  the  universal  wavefunction  Y.  In  the  MIW 
approach,  such  a  branch  will  correspond  to  a  subset  Xe  C  X  of  world-particles,  still  extremely  large 
in  number,  that  are  close  together  (on  a  macroscopic  scale)  in  configuration  space  and  mostly  far  (on 
a  macroscopic  scale)  from  the  rest  of  X.  Thus,  to  a  good  approximation,  this  set  of  world-particles  will 
evolve  autonomously,  on  the  time  scale  of  interest. 

Now  consider  a  particular  world-particle,  xj  e  Xp,  and  a  particular  part  of  the  vector  xy, 
comprising  a  lower-dimensional  vector  by,  that  contains  the  variables  that  encode  the  macroscopic 
fact  of  the  decision  by  an  experimenter.  Bob,  about  what  experiment  to  perform.  That  is,  the  decision 
is  the  same  in  every  world  under  consideration  here.  However,  there  is  no  reason  in  the  theory  for 
by  to  be  correlated  with  the  other  variables  in  xy,  until  the  experiment  is  actually  performed.  Once  it 
is  performed,  the  value  of  by  will  have  a  direct  (second-order  in  time)  effect  on  the  vicinal  (in  the  3D 
sense)  variables  in  xy,  say  By.  Then,  in  the  MIW  approach,  the  change  in  By  will  cause  a  change  (again 
second-order  in  time)  in  xj-  G  Ny  (note  that  Ny  is  only  a  tiny  subset  of  Xe).  Now  say  that  in  Ye,  there  is 
entanglement  between  the  system  corresponding  to  By,  and  a  distant  system  corresponding  to  some 
other  observables  Ay.  Then  the  change  in  x/;.  will  be  a  change  not  just  in  B/.  but  also  in  A/,.  This  is 
because  the  MIW  acceleration  function  a^iw  pays  no  heed  to  vicinality  in  the  3D  sense — it  cares  only 
about  the  fact  that  the  worlds  x/.  G  Ny  are  in  the  vicinity  of  the  world  xy  in  configuration  space. 
Finally,  the  same  interaction  will  cause  a  (second-order  in  time)  change  in  Ay,  the  distant  observables 
in  the  original  world  j,  since  for  at  least  some  of  the  x/c  E  Ny  we  will  have  xy  E  N/c. 
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Thus  it  is  clear  in  principle  how  the  MIW  approach  can  give  rise  to  nonlocality  in  Bell's  1964 
sense,  with  Bob's  choice  bj  causing  a  change  in  the  variables  Ay,  variables  which  correspond  to  the 
outcome  observed  by  Alice  in  her  distant,  perhaps  even  space-like  separated,  lab,  all  in  the  same 
world  j.  It  might  seem  that  Bob  could  use  this  to  signal  to  Alice  faster  than  light.  However,  that  is  not 
so,  for  essentially  the  same  reason  that  signal-locality  is  respected  in  Bohmian  mechanics.  Alice  and 
Bob  are  ignorant  of  which  world  j  they  inhabit.  They  know  only  the  macroscopic  configuration  Xf  in 
which  their  world  must  lie.  A  successful  MIW  theory  would  thus  exhibit  nonlocality  at  the  hidden 
level  of  an  individual  world,  but  not  at  the  observable  level,  averaging  over  all  the  worlds  in  Xp. 

4.2.  Simulations 

Can  we  explicitly  simulate  this  in  our  MIW  approach,  showing,  for  example,  the  violation  of 
a  Bell  inequality,  or  even  just  the  EPR  paradox,  without  signalling?  Here  we  run  up  against  the 
limitations  in  the  development  of  our  approach. 

In  our  2014  paper  [29],  we  introduced  a  toy  MIW  theory  to  describe  a  very  simple  universe, 
comprising  a  single  particle  in  ID.  We  showed  analytically  that,  in  the  limit  of  a  large  number 
of  worlds,  it  gave  the  correct  ground  state  distribution  for  a  harmonic  oscillator,  and  the  correct 
description  for  the  first  and  second  moments  of  free  particle  evolution.  We  showed  numerically  that 
it  could  reproduce,  at  least  qualitatively,  the  double-slit  interference  phenomenon.  Furthermore,  we 
argued  that  it  can  plausibly  reproduce  other  generic  quantum  phenomena  such  as  barrier  tunnelling 
and  reflection. 

To  demonstrate  nonlocal  correlations  between  measurements  on  distant  quantum  systems,  we 
obviously  require  more  than  one  particle.  One  might  think  one  would  require  particles  to  model 
the  measuring  apparatuses  as  well  as  having  at  least  two  entangled  particles.  However,  since  the 
position  of  a  particle  is  reified  in  the  MIW  approach,  similarly  to  Bohm's  theory  [4,5],  a  particle  can 
always  represent  its  own  measured  position  (this  differs  from  its  momentum,  as  any  measurement 
of  momentum  will  rely  on  the  action  of  the  quantum  potential  in  Bohm's  theory,  or  the  interworld 
interaction  potential  in  our  approach).  Moreover,  if  we  consider  a  particle  on  a  spring,  with  an 
externally  controllable  spring  constant  k,  then  either  its  initial  position  or  its  initial  momentum 
can  be  encoded  in  position  at  a  nominated  time.  This  is  because  such  evolution  turns  momentum 
information  into  position  information  at  time  T/4  =  (tc/2)  \/'m / k,  and  back  into  position  information 
at  time  T  /  2.  (Alternatively,  the  same  can  be  achieved  even  by  free  evolution  with  an  externally 
controlled  mass  m,  an  assumption  which  does  not  violate  any  of  the  MIW  framework.  By  reducing  the 
mass,  the  time  taken  for  the  momentum  to  be  encoded  into  the  position  can  be  made  arbitrarily  small; 
by  increasing  the  mass,  the  position  remains  as  it  was  initially  to  an  arbitrarily  good  approximation. 
Thus,  either  position  or  momentum  could  be  measured  at  a  nominated  time.)  In  this  way,  nonlocal 
correlations  can  be  explored  with  a  world  of  just  two  ID  particles.  This  scenario  could  be  used  not 
only  to  demonstrate  the  EPR  paradox,  but  also  Bell  nonlocality,  as  proven  by  Bell  himself  in  1986  [37]. 

Recently,  we  have  generalised  MIW  for  a  2D  particle,  which  is  equivalent  to  two  ID 
particles,  with  some  numerical  success  [38].  However,  that  study  was  restricted  to  finding  ground 
states — stationary  states  with  no  nodes  in  the  wavefunction.  The  same  paper  did  consider  finding 
the  first  excited  states  for  a  ID  particle,  but  only  for  symmetric  potentials  for  which  the  node  at  x  =  0 
can  be  put  in  by  hand.  This  restriction  was  necessary  because,  it  is  now  clear,  the  toy  MIW  model  of 
Reference  [29]  cannot  reproduce  stationary  excited  states  at  all.  More  generally,  it  cannot  be  expected 
to  quantitatively  reproduce  dynamical  quantum  evolution  in  which  nodes  appear  and  disappear.  In 
the  proposal  of  Bell  [37],  to  violate  a  Bell  inequality  by  free  evolution  of  two  entangled  particles,  the 
required  initial  wavefunction  has  two  nodes  and  can  be  expected  to  develop  more  during  evolution. 
Thus,  to  study  nonlocality  in  the  MIW  approach  we  certainly  need  to  go  beyond  the  toy  model  of 
Reference  [29].  Progress  in  this  direction  is  the  topic  of  the  next  two  sections. 
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5.  MIW  Beyond  the  Toy  Model 

The  MIW  toy  model,  for  a  nonrelativistic  ID  quantum  particle  of  mass  m  moving  in  a  potential 
V(q),  has  N  worlds  each  containing  a  ID  particle.  Denoting  the  positions  and  momenta  of  these 
world-particles  by  X  =  (xy,  X2, . . . ,  Xjy)  and  P  =  (pi,  pi>  ■  ■  ■  ,  P N ) ,  their  motion  is  described  in  the  toy 
model  by  a  Hamiltonian  of  the  form  [29] 

H(X,P)  =  £  +  V(xn)j  +  Utoy(X).  (3) 


The  first  term  is  just  the  sum  of  N  classical  Hamiltonians,  one  for  each  world.  It  is  the 
second  term,  a  potential  energy  responsible  for  interactions  between  the  worlds,  which  accounts  for 
quantum  phenomena: 


Iftoy(X) 


h  2_  *  /  1 _ 1  \ 

8m  „=1  V  xn+ 1  -  Xn  Xn  -  Xn- 1 ) 


(4) 


For  evaluation  purposes  we  formally  define  Xo  —  oo  and  xpj+i  =  oo,  and  the  ordering 

xi  <  x2  <  ■  ■  ■  <  x,v  has  been  assumed  (this  ordering  is  preserved  under  time-evolution,  due  to 
the  repulsive  nature  of  the  potential).  The  motion  of  each  world-particle  is  determined  by  the  usual 
Hamiltonian  equations  of  motion,  i.e.. 


dH  _i 
Xn  =  5 —  =  rn  yn/ 
dpn 


Pn 


rl !  i 

dxn 


-v'{xn) 


r)  ft  toy 

dxn 


(5) 


Note  that  the  interworld  potential  Ut0y(X)  vanishes  in  the  classical  limit  h  —  0,  and  also  for 
the  case  N  —  1.  In  either  of  these  cases,  each  world  evolves  independently,  according  to  Newton's 
laws.  More  generally,  however,  for  h  ^  0  and  N  >  1  the  interworld  potential  lft0y(X)  leads  to  forces 
on  each  world  that  act  to  reproduce  quantum  phenomena  such  as  Ehrenfest's  theorem,  spreading 
of  wave  packets,  tunneling  through  a  barrier,  and  interference  effects  [29].  For  the  case  of  a  ID 
oscillator,  V (q)  —  \moj2q2,  it  has  further  been  shown,  in  the  limit  N  —?  oo,  that  the  average  energy 
per  world  of  the  MIW  ground  state  converges  to  the  quantum  groundstate  energy  jtuo  [29],  and  that 
the  corresponding  stationary  distribution  of  worlds  samples  the  usual  quantum  Gaussian  probability 
distribution  [29,39]. 

The  form  of  Lftoy(X)  in  Equation  (4)  above  is  a  sum  of  three-body  terms,  leading  to  a  force 
on  the  w-th  world  that  depends  on  the  positions  of  the  two  neighbouring  worlds  on  either  side. 
However,  while  this  form  is  sufficient  to  reproduce  the  quantum  phenomena  noted  above,  we  have 
found  that  it  is  too  simple  to  model  the  behaviour  of  quantum  wave  functions  with  nodes.  Hence, 
as  noted  in  the  previous  section,  we  must  turn  to  more  complex  forms  of  the  interworld  potential, 
involving  interaction  between  greater  numbers  of  neighbouring  worlds. 

Fortunately,  there  is  a  great  deal  of  freedom  in  choosing  this  potential  in  the  MIW  approach. 
It  is  possible  this  freedom  could  be  curtailed  via  suitable  physically-motivated  axioms  (such  as,  for 
example,  requiring  Ehrenfest's  theorem  to  hold).  Here,  however,  we  take  a  nonaxiomatic  approach, 
to  show  how  the  interworld  potential  of  the  toy  model  can  be  straightforwardly  generalised  to 
allow  direct  interactions  between  an  arbitrary  number  of  worlds,  in  a  manner  corresponding  to 
greater  accuracy  in  approximation  of  the  Bohmian  acceleration.  This  greater  accuracy  supports  a 
corresponding  expectation  of  being  able  to  successfully  model  wave  function  nodes. 

In  particular,  it  was  shown  in  section  II. D  of  Reference  [29]  that  a  suitable  general  form  for  the 
interworld  potential  in  the  ID  case  is  obtained  by  replacing  Ut0y(  X)  in  the  Hamiltonian  (3)  by 


U(X) 


(6) 
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where  P„  and  P'  are  approximations  to  a  probability  density  P(cj)  at  q  —  xn,  and  where  P(q) 
corresponds  to  some  smoothing  of  the  empirical  distribution  of  the  world  positions  { X\ ,  x^,  ■  ■ .  ,Xjy}. 
Equation  (6)  may  be  regarded  as  an  approximation  of  the  quantum  potential  in  Bohm's  theory  [4,5] 
(for  which  P(q )  =  Y (q)  2),  and  hence  is  expected  to  reproduce  quantum  evolution  more  and  more 
closely  as  this  approximation  is  improved.  The  toy  model  potential  in  Equation  (4)  is  the  simplest 
such  approximation  [29]. 

To  go  to  higher  order  approximations,  we  tried  two  different  methods.  The  first  we  call  the 
rational  smoothing  method  and  the  second  the  equivariance  method.  In  rational  smoothing  we  take 
a  systematic  approach  to  approximate  P(q)  and  its  derivative  by  ratios  of  polynomials,  where  the 
order  of  these  polynomials  determines  the  number  of  neighbours  that  each  world  directly  interacts 
with.  This  approach  is  developed  in  the  next  subsection,  followed  by  an  example.  It  will  be  applied 
to  the  description  of  nodes  in  Section  6.2.  rational  smoothing  is  less  computational  resource-intensive 
compared  to  the  equivariance  Method.  The  latter  is  discussed  in  Appendix  A. 


5.1.  Constructing  Generalised  Interworld  Potentials 

If  P(q)  is  some  smooth  probability  density  as  above,  which  approximates  the  distribution  of 
world  positions,  {x\,xi, . . . ,  Xjy },  then  the  cumulative  distribution  C(x)  f  dq  P(q)  must  satisfy 
(at  least  approximately) 

C(xn)  =  dq  P(q)  =  un  :=  (7) 

Thus,  the  area  under  P(q)  is  divided  into  N  neighbourhoods,  each  containing  one  world  and 
having  area  1/N.  It  will  be  assumed  for  simplicity  that  P(q)  is  nonzero  almost  everywhere,  implying 
that  C(x)  is  strictly  monotonic  and  hence  invertible. 

To  systematically  approximate  the  quantum  force  at  xn  to  a  given  accuracy,  we  first  define  the 
inverse  of  the  cumulative  distribution  by  y(u)  C  ‘  1  (u).  It  follows  immediately  from  Equation  (7) 
that 

Xn  —  y{Un).  (8) 

Further,  differentiating  u  =  C(y)  with  respect  to  y  and  using  d/dy  =  (1  / y')(d/du),  gives  P(y)  = 
1  /y',  P'(y)  —  —  y"/(y')3,  and  hence  that 


P'(y )  _  y" 
P(y)  (yT 


(9) 


Our  aim  is  now  to  approximate  this  last  expression  by  a  function  of  the  world  positions,  and 
then  substitute  this  approximation  into  the  right  hand  side  of  Equation  (6)  to  obtain  a  corresponding 
form  for  the  interworld  potential  LI (X). 

In  particular,  expanding  y(un+c)  in  a  Taylor  series  about  y{un)  gives  the  approximation 


xn+c  xn  —  y(un+c)  1 l{un)  ~  YL  J|  (]\[)  ^ 


(10) 


to  accuracy  0(1  /NL),  where  y„  denotes  the  Zth  derivative  of  y ( u )  at  u  —  un ■  The  first  L  derivatives 
of  y  can  therefore  be  approximated  to  this  accuracy  by  choosing  a  set  of  coefficients  { ac/ }  such  that 


Y^otdC1'  =l\6lv,  1,1'  =  1,2,...,  L.  (11) 


It  is  shown  how  to  construct  suitable  { ac/ }  below.  Equations  (10)  and  (11)  immediately  yield 
the  approximations 


t/(,) 

N1 


E«c/(*n+c  xn)r 


(12) 
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accurate  to  0(1/ NL),  for  l  —  1,2 Finally,  substitution  into  Equation  (9)  and  hence  into 
Equation  (6)  gives  the  corresponding  interworld  potential. 


N 


U*,L(X)  =  E  U«,L,n(X) 

n—1 


ft  r— '  j  Lc<Xc2(Xn+c  Xu)  1 
8m„=lUEcM*n+c-*n)]2J 


(13) 


Note  that  this  interworld  potential  is  translation  invariant,  and  scales  as  1  /  A2  under  x„  — »■ 
Axn.  It  follows  that  analogues  of  Ehrenfest's  theorem  and  wavepacket  spreading  hold  for  all  such 
potentials  [29]. 

An  advantage  of  Hamiltonian  based  methods  is  that  we  can  use  a  symplectic  numerical 
integrator  to  make  sure  that  the  total  energy  is  conserved  [40].  The  quantum  force  for  the  mth 
world-particle  can  be  calculated  from  Equations  (5)  and  (13)  as: 


dUgjJXj  =  _  £  31WX) 

dxm  n  | 


(14) 


5.2.  Examples 

To  obtain  an  explicit  example  of  the  interworld  potential  in  Equation  (13),  let  M  denote  an  L  x  C 
matrix  with  coefficients  M/c  =  c1 ,  and  A  denote  the  C  x  L  matrix  with  coefficients  xcj,  where  c  ranges 
over  some  set  of  C  integers.  Equation  (11)  can  then  be  written  in  the  matrix  form 

MA  =  A  :=diag[l!,2!,...,L!|.  (15) 


The  existence  of  a  solution  requires  C  >  L.  Further,  since  the  values  of  dc\  and  otc 2  are  required 
in  Equation  (13),  one  must  have  L  >  2.  The  corresponding  solution  is  then  A  =  M_1A,  where  M”1 
denotes  the  inverse  of  M  for  C  —  L,  and  a  pseudo-inverse  of  M  for  C  >  L. 

It  follows  that  the  simplest  interworld  potential  constructed  in  this  way  corresponds  to  C  —  L  = 
2.  Labelling  the  two  values  of  c  by  c  —  ±1  corresponds,  via  Equation  (12),  to  approximating  the 
derivatives  of  y{xn)  via  the  values  of  the  nearest  neighbours  xn±i-  Ordering  the  values  of  c  as  —1, 1, 
the  corresponding  coefficients  { txci }  in  Equation  (13)  are  then  given  by 


A  =  M”1  A  = 


1  0 
0  2 


-1/2  1  \ 
1/2  1  J  ‘ 


(16) 


Surprisingly,  this  is  not  actually  equivalent  to  the  toy  model  of  Reference  [29],  even  though  it 
involves  the  same  number  of  neighbours  in  the  potential. 

The  simplest  higher-order  interworld  potential  corresponds  to  C  =  L  =  3,  but  odd  values  of  C 
necessarily  introduce  an  unphysical  left-right  asymmetry,  with  each  world  being  coupled  to  different 
numbers  of  neighbouring  worlds  on  either  side  via  l lK,L,n{X).  To  preserve  symmetry  we  therefore 
next  consider  the  case  C  —  L  —  4.  Labelling  the  four  values  of  c  by  c  —  ±1,  ±2,  the  derivatives  of 
y{xn)  in  Equation  (12)  can  be  calculated  with  values  of  the  nearest  and  next  nearest  neighbours  xn±\ 
and  xn±2-  If  we  order  the  values  of  c  as  —2,  —1, 1,2,  the  coefficients  { txcj }  in  Equation  (13)  are  then 
given  by 
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(17) 


6.  Numerical  Results 

To  test  our  higher-order  methods  against  the  original  toy  model,  we  apply  them  to  the  problem 
of  finding  the  ground  and  the  first  excited  state  of  a  harmonic  oscillator.  It  is  the  latter,  containing  a 
node,  where  the  toy  model  fails  and  it  is  necessary  to  use  a  higher-order  interworld  potential.  In  this 
section  we  use  dimensionless  configuration  coordinates 


Xn  :=  V2mca/h  xn 

(18) 

and  dimensionless  times 

T  oot/ln, 

where  co  is  the  harmonic  oscillator  frequency. 

6.1.  Toy  Model 

(19) 

Figure  lb  shows  the  result  of  applying  the  toy  model  in  Equations  (3)  and  (4)  to  the  ground 
state  of  a  harmonic  oscillator,  corresponding  to  V (q)  —  jot (02x2.  In  this  test  we  distributed  50 
world-particles  with  xn  determined  by  inverting  Equation  (7)  for  the  groundstate  probability  density 
Pf  °  (q)  (Figure  la).  We  then  evolved  these  under  the  single-particle  harmonic  oscillator  potential 
\mcv2x2  and  the  quantum  interworld  potential.  We  see  that  as  expected  for  a  stationary  quantum 
potential,  the  classical  and  quantum  forces  cancel  each  other  and  the  world-particles  stay  stationary. 
The  slight  oscillations  for  the  world-particles  near  the  boundary  is  due  to  differences  between  the 
toy-model  and  Bohmian  potentials  in  areas  with  high  curvature  and  low  sampling.  These  differences 
imply  an  exact  stationary  state  for  the  MIW  potential  that  has  slightly  different  values  for  the  xn  for 
any  finite  value  of  N  [29]. 

However,  if  we  distribute  world-particles  based  on  the  probability  density  of  the  first  excited 
state  of  a  harmonic  oscillator  (Figure  2a)  and  apply  the  same  model,  the  world-particles  will  not  stay 
stationary  as  expected  (Figure  2b).  We  also  implemented  the  nearest-neighbour  interworld  potential 
defined  by  Equation  (17),  and  found  very  similar  behaviour  to  the  original  toy  model.  That  is,  it  fails 
to  support  stationary  configurations  corresponding  to  excited  energy  eigenstates. 
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Figure  1.  (a)  50  world-particles  (that  is,  50  different  worlds  describing  a  single  particle)  are  distributed 
based  on  the  probability  density  of  the  ground  state  of  a  harmonic  oscillator,  (b)  Trajectories  of  the 
50  world-particles.  As  expected,  the  classical  and  quantum  forces  approximately  cancel  each  other 
and  the  world-particles  stay  approximately  stationary.  The  slight  oscillations  near  the  boundary  are 
due  to  approximation  of  the  Bohmian  potential  by  the  toy-model  potential. 


The  unsuccessful  simulation  results  in  Figure  3  show  that  the  problem  of  the  nodes  cannot  be 
fixed  by  increasing  the  sampling.  In  the  case  of  a  node,  because  the  probability  density  is  zero,  there  is 
always  a  second  order  curve  between  the  two  world-particles  adjacent  to  the  node  that  can  never  be 
correctly  estimated  by  toy  model  approximation.  This  leads  to  a  poor  approximation  of  the  Bohmian 
force  for  the  world-particles  near  the  node  and  instead  of  staying  stationary  the  world-particles  move 
towards  the  gap  and  fill  the  gap. 

The  general  problem  with  areas  of  low  probability  density,  and  in  particular  in  the  region  of 
the  node  of  the  quantum  state  at  x  —  0  in  Figure  2a,  is  that  sampling  is  low  and  nearest-neighbour 
approximations  are  not  valid  for  calculating  the  quantum  potential.  If  the  probability  is  low  but 
non-zero,  theoretically,  we  can  increase  the  number  of  worlds  until  we  reach  a  good  sampling  in  those 
areas.  To  test  the  simulation  with  higher  sampling,  we  tried  5000  worlds.  To  reduce  the  simulation 
time,  we  only  focused  on  the  20  world-particles  around  the  node  (Figure  3a).  To  apply  the  correct 
boundary  condition  for  truncated  area,  we  kept  five  world-particles  near  each  boundary  artificially 
fixed. 


Figure  2.  (a)  40  world-particles  are  distributed  based  on  the  probability  density  of  the  first  excited 
state  of  a  harmonic  oscillator,  (a)  Trajectories  for  the  first  excited  state  of  a  harmonic  oscillator  using 
the  potential  in  Equation  (4).  The  world-particles  do  not  stay  stationary.  Particularly  those  near  the 
node  move  towards  the  middle  and  fill  the  gap. 
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Figure  3.  (a)  5000  world-particles  are  distributed  based  on  the  probability  density  of  the  first  excited 
state  of  a  harmonic  oscillator.  Only  20  world-particles  around  the  node  in  the  middle  are  shown, 
(b)  Trajectories  of  the  world-particles  for  the  initial  distribution  in  (a).  To  apply  the  correct  boundary 
condition,  we  kept  five  world-particles,  near  each  boundary,  fixed  and  simulated  the  dynamics  of  the 
remaining  10  world-particles  in  the  middle,  time-steps  are  10~8.  Since  the  nearest  neighbours  of  the 
node  do  not  stay  near  the  starting  point  and  move  to  the  middle  of  the  node  after  approximately  0.1  of 
a  period,  we  did  not  continue  the  simulation  for  a  full  period,  (c)  The  same  test  as  (b)  with  time-steps 
of  10~9.  Thus,  the  failure  of  the  simulation  is  not  an  artefact  of  large  time-steps. 


6.2.  Higher  Order  Potential 

To  test  a  higher  order  approximation,  we  applied  the  rational  smoothing  model  in  Equation  (13) 
for  the  case  L  =  4  (equivalent  to  a  5-world  interaction),  to  the  first  excited  state  of  a  harmonic  oscillator 
with  5000  worlds.  This  corresponds  to  the  example  in  Equation  (17). 

In  the  first  test  we  kept  all  the  worlds  artificially  stationary  except  for  the  two  worlds  around  the 
node.  Figure  4  shows  that  these  two  middle  worlds  stay  stationary  which  means  that  the  quantum 
potential  for  those  is  well  approximated. 


Postion  X 


Figure  4.  Simulation  of  only  two  worlds  adjacent  to  the  node  for  the  first  excited  state  of  a  harmonic 
oscillator.  The  initial  positions  were  set  by  considering  5000  worlds  in  total  to  describe  the  excited 
state.  To  apply  the  boundary  condition,  the  rest  of  the  world-particles  were  kept  fixed.  For  the 
evolution,  rational  smoothing  with  L  =  4  (equivalent  to  5-world  approximation)  is  used.  Time  step  is 
10 ~9.  The  two  world-particles  stay  stationary  as  expected. 
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In  the  second  test  we  tried  the  same  scenario  but  this  time  simulated  10  worlds  around  the  node 
and  kept  the  rest  stationary.  Figure  5  shows  that  the  simulation  was  successful  in  the  sense  that 
world-particles  stayed  close  to  the  starting  position  and  did  not  move  towards  the  gap  i.e.,  the  low 
density  region,  corresponding  to  a  node  of  the  associated  excited  state,  is  preserved  by  the  interaction, 
and  the  evolution  is  approximately  stationary.  The  obvious  difference  compared  to  Figure  4  is  the 
oscillations  around  the  starting  point.  The  oscillations  appear  chaotic,  as  might  be  expected  for  a 
system  of  nonlinearly  coupled  harmonic  oscillators.  These  oscillations  might  also  be  due  to  the  time 
step  being  too  long.  We  were  not  able  to  test  this  conjecture  because  the  test  in  Figure  5  took  a  few 
days  on  our  desktop  computer  and,  due  to  time  constraints,  we  could  not  run  it  with  smaller  time 
steps. 

We  repeated  the  same  test  with  7-worlds  approximation  (L  =  6).  Figure  6  shows  that  this  higher 
order  approximation  decreases  oscillation  compared  to  the  results  in  Figure  5.  Hence,  the  simulation 
of  a  wave  function  with  nodes,  which  failed  for  the  toy  model  (Figures  2  and  3),  is  seen  to  become 
more  and  more  accurately  modelled  in  the  MIW  approach  as  the  number  of  directly  interacting 
worlds  is  increased  (Figures  5  and  6). 

Similar  convergence  might  also  be  possible  with  the  equivariance  Method  presented  in 
Appendix  A.  However,  the  results  of  the  simplest  test,  with  a  potential  involving  five  interacting 
worlds,  were  not  as  positive  as  those  of  the  rational  smoothing  explored  in  this  section.  They  are  also 
more  numerically  intensive,  so  we  did  not  pursue  it  further. 


Figure  5.  Simulation  of  10  worlds  (rather  than  two,  as  in  Figure  4),  five  on  either  side  of  the  node  for 
the  first  excited  state  of  a  harmonic  oscillator.  Other  details  are  as  in  Figure  4. 
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Figure  6.  Simulation  of  the  first  excited  state  of  a  harmonic  oscillator  with  rational  smoothing,  L  =  6, 
(7-world  approximation).  Other  details  are  as  in  Figure  5.  The  oscillations  are  much  smaller  than 
those  in  the  L  =  4  approximation  in  Figure  5. 

7.  Conclusions 

Bell's  theorem  of  1964  showed  that  any  deterministic  interpretation  of  quantum  mechanics  must 
be  nonlocal  (Section  2).  Bohm's  theory  of  1952  is  the  example  par  excellence  of  a  nonlocal,  deterministic 
theory  (Section  3).  The  theoretical  approach  we  introduced  in  2014 — Many  Interacting  Worlds — is 
also  deterministic  and,  if  it  is  to  succeed  in  replicating  Bohmian  mechanics  (and  thus  all  quantum 
phenomena),  must  be  nonlocal  (Section  4).  It  is  conceptually  clear  how  an  interworld  potential  can 
lead  to  nonlocality  in  Bell's  1964  sense,  and,  perhaps  surprisingly,  this  could  in  principle  be  simulated 
with  a  universe  of  only  two  particles,  in  ID  (one  spatial  dimension),  if  an  externally  controlled 
spatially  localized  potential  is  allowed,  as  discussed  in  Section  4. 

Unfortunately,  this  proof-of-principle  simulation  of  Bell-nonlocality  cannot  be  done  using  the  toy 
model  for  an  interworld  potential  introduced  in  Reference  [29].  The  reason  is  that  it  cannot  deal  with 
wavefunctions  having  nodes,  and  nodes  are  certainly  necessary  for  modelling  the  entangled  states 
and  measurements  necessary  for  violating  a  Bell  inequality  [37].  When  one  prepares  a  distribution 
of  worlds,  in  one  dimension,  corresponding  to  a  stable  excited  quantum  state,  the  dynamics  of 
the  toy  model  causes  the  gap  between  worlds  where  the  node  should  collapse,  and  no  stationary 
configuration  is  reached.  Dealing  with  nodes  is  a  problem  in  many  quantum  simulation  methods 
based  on  Bohmian  mechanics  [41].  Nodes  should  not  be  a  problem  for  interpretations  involving  a 
continuum  of  worlds  [32-35],  as  they  are  formulated  to  be  exactly  equivalent  to  quantum  mechanics. 
However,  as  remarked  in  Section  3,  our  view  is  that  these  interpretations  do  not  solve  the  conceptual 
problems  of  the  Everettian  many-worlds  interpretation. 

Here  we  showed  that  this  problem  of  nodes  in  our  discrete  MIW  approach  may  not  be 
fundamental,  but  rather  may  be  a  result  of  using  a  too  simple  form  for  the  inter-world  interaction 
potential.  By  using  a  higher-order  approximation  to  define  our  interworld  potential  from  Bohm's 
quantum  potential,  we  were  able  to  show  that  a  gap  in  the  configuration  of  worlds,  corresponding 
to  the  node  of  the  first  excited  harmonic  oscillator  energy  eigenstate,  can  remain  open  for  at  least  a 
full  harmonic  oscillator  period,  and  perhaps  indefinitely.  The  world  configurations  were  not  exactly 
stationary,  but  rather  had  high  frequency  irregular  oscillations.  However,  by  increasing  the  order  of 
the  approximation,  the  size  of  the  oscillations  could  be  reduced. 
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Our  simulations  considered  only  the  dynamics  near  the  node,  with  more  distant  worlds 
artificially  held  fixed.  Whether  the  node  would  remain  stable  if  all  worlds  were  allowed  to  evolve 
according  to  the  MIW  dynamics  is  an  open  question.  In  addition,  our  simulations  were  restricted  to 
one  particle  in  ID.  The  MIW  approach  has  been  successfully  used  to  simulate  one  particle  in  2D  (or, 
equivalently,  two  particles  in  ID)  [38,42]  and  3D  [42],  but  only  to  find  the  ground  state  configuration. 
Combining  these  research  directions  to  be  able  to  simulate  stable  excited  states  for  two  particles  in 
ID,  and  beyond,  is  a  challenge  for  future  work.  Finally,  to  realise  a  simulation  of  the  Bell  experiment 
described  in  Section  4  would  require  simulating  not  just  stationary  nodes  but  also  dynamical  nodes, 
that  may  appear  and  then  disappear  in  an  instant,  which  are  notoriously  difficult  to  deal  with  in 
Bohmian-inspired  numerical  approaches  [41].  Thus,  much  work  remains  to  be  done,  but  the  positive 
results  reported  here  are  encouraging. 
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Abbreviations 

The  following  abbreviations  are  used  in  this  manuscript: 

MIW  Many  Interacting  Worlds 

EPR  Einstein,  Podolsky,  and  Rosen 

ID  One  (spatial)  dimension 

2D  Two  (spatial)  dimensions 

3D  Three  (spatial)  dimensions 

Appendix  A.  A  Equivariance  Method 

Here  we  introduce  an  alternative  higher  order  method  to  rational  smoothing,  motivated  by  the 
equivariance  property  of  Bohmian  mechanics  [11].  In  particular,  for  a  given  configuration  of  worlds, 
{  x  ] ,  i'2, . . . ,xn) ,  we  construct  a  smooth  polynomial  probability  density,  Pn  (x),  in  the  region  of  each 
world  xn,  and  use  this  to  approximate  the  Bohmian  potential.  The  coefficients  of  the  polynomial 
are  determined  by  requiring  that  area  under  Pn(x )  between  xn  and  xn+x  is  equal  to  the  constant 
value  1/N  for  each  world,  analogously  to  Equation  (7).  The  accuracy  of  this  equal-probability  or 
"equivariance"  method  will  increase  with  the  degree  of  the  polynomial. 

Here  we  will  illustrate  the  equivariance  method  for  the  case  where  P„  (x)  is  third-order, 
corresponding  to  direct  interactions  between  sets  of  five  adjacent  worlds.  Thus, 

P„  (x)  =  an  +bnx  +  cnx2  +  dnx3.  (20) 

The  equivariance  method  then  requires  that 

fxn+l  1 

/  Pn(x)dx  —  const.  =  N  .  (21) 

J  Xn 

To  determine  the  coefficients  of  P„  (x)  in  Equation  (20),  we  use  four  equations  below  based  on 
the  positions  of  the  five  worlds,  xn-2,  x„_i,  xn,  xll+\  and  xn+2' 


fX”  1  Pn(x)dx  —  N-1  (22) 

J  Xn—2 

[ X"  Pn(x)dx  —  AT1  (23) 

Jxn-\ 

fxn+ 1 

/  Pn(x)dx  —  N-1  (24) 

Jxn 

rxn+  2 

/  Pn{x)dx  =  N~ 1  (25) 

Jxn+l 


Substituting  P„  (x)  from  Equation  (20)  and  evaluating  the  integrals  we  get 

an  (*„_!  -  x„_2)  +  bn  -  xl-l)  +  Cn  UXl- 1  “  xl-l)  +  l(Xi-l  ~  Xn-l)  =  N _1 

a„  (xn  —  xn_i)  +  bn  l(x 2  -X2_x)  +  c„  l(x3  —  x3_x)  +  dn \{x\  -  x\_x)  = 
an  (xn+i  -xn)+bn  l(xl+1  ~xl)+cn  l(x3n+l  -x3)  +dn  \(xfl+1  -  x\)  = 

«»  (Xn+ 2  ~  *H+l)  +  bn  2(4+2  ~  4+1 )  +  3  (4+2  ~  4+1 )  +  dn  1(4+2  ~  4+1 )  =  N”1' 


(26) 
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If  we  define  matrix  Kn  as 


K„  = 


(  0«_i  -  xn-2) 

( xn  xn—p) 
(xn+ 1  —  ^n) 
\  ixn+ 2  —  *n+l) 


we  can  rewrite  Equation  (26)  as 


44-1  “4-2)  44-1  ~4-2) 
44“4-l)  44“  4-1 ) 


44+i~4) 


2  V  n+1 

44+2“  4+1 ) 


3(*n+l-*n) 

3  (4+2  ~  4+ 1 ) 


44-1  ~4-2)  \ 
44-4-4 
1(4+1  -4) 
1(4+2- 4+i )  / 


(27) 


Kn 


/  tin  ^ 

{ N_1  \ 

bn 

N-1 

Cn 

AT1 

\  dn  j 

l  N ■ 1  / 

(28) 


Therefore: 


(  \ 


un 

Cn 


=  K~' 


(  N~1  \ 
AT1 

TV  r— 1 


(29) 


We  used  Matlab  to  invert  Kn  and  find  coefficients  an,bn,cn,  and  dn.  The  resulting  equations 
were  too  lengthy  to  include  here.  Substituting  these  coefficients  into  Equation  (20),  we  can  find 
the  probability. 

To  evaluate  the  quantum  force,  we  use  an  interworld  interaction  potential  of  the  form  of 
Equation  (6)  and  Pn  as  in  Equation  (20)  (with  an,  bn,  c„  and  dn  taken  from  (29))  to  give 

jj h  r— i  (  bn  +  2cnx„  +  3dnx^  \  (30) 

8tn  \an  +  bnxn  +  CnXpj  +  dnxfi ) 

This  may  be  compared  with  the  corresponding  interaction  potential  U^i  in  Equation  (13) 
obtained  via  rational  smoothing.  The  quantum  force  on  each  particle  is  evaluated  similarly  to 
Equation  (14).  We  used  Matlab  to  evaluate  the  analytical  derivatives  of  these  terms  for  the  simulation. 
The  resulting  equation  is  too  lengthy  to  include  here  (c.  360,000  characters  for  N  —  5). 

Figure  7  shows  the  result  of  applying  the  equivariance  method  to  the  first  neighbours  of  the  node 
in  the  first  excited  state  of  a  harmonic  oscillator.  It  shows  that  these  world-particles  do  not  move  into 
the  gap  from  their  original  position,  but  rather  undergo  some  oscillatory  motion.  This  is  poorer 
behaviour  than  the  corresponding  rational  smoothing  simulation  shown  in  Figure  4.  The  simulations 
were  also  slower  because  of  the  complexity  of  the  analytical  form  of  the  force  law,  mentioned  above. 
For  these  reasons  we  have  not  pursued  this  method  further. 


Version  September  6,  2018  submitted  to  Entropy 


18  of  19 


position,  x 


Figure  7.  The  first  excited  state  of  a  harmonic  oscillator  is  simulated  using  5-world  approximation  in 
equivariance  Mmethod.  Five-thousand  worlds  are  used  and  only  the  two  world-particles  next  to  the 
node  are  simulated.  The  rest  are  kept  stationary,  similar  to  Figure  4  for  the  rational  smoothing  case. 
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